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SERVICES  GROUP 

1  July  1970 


Headquarters 

United  States  Air  Force 

AFTAC/VSC 

Washington,  D.  C.  20333 
Attention:  Lt.  John  Woods 

Subject:  Third  Quarterly  Report  Covering  Period  1  March  1970 

to  31  May  1970 

Identification:  AFTAC  Project  No.  :  VELA/T/0702/B/ASD 
Project  Title:  Statistical  Discrimination 
ARPA  Order  No. :  624 
ARPA  Program  Code:  9F10 

Name  of  Contractor:  Texas  Instruments  Incorporated 
Contract  Number:  F33657-70-C-0311 
Effective  Date  of  Contract:  25  August  1969 
Amount  of  Contract:  $222,  250 
Contract  Expiration  Date:  14  September  1970 
Project  Manager:  Stanley  J.  Laster 


Gentlemen: 

MAJOR  ACCOMPLISHMENTS  (TASKS  REFER  TO  WORK  STATEMENT) 

(1)  ^hort-period  event  classification  studies  are  essentially 
complete  except  for  the  processing  of  thirty  selected 
events  at  SAAC  and  eight  events  supplied  by  the  contract 
monitor. The  Seismic  Event  Classification  Software 
Package  has  been  modifed  to  accept  long -period  statistics. 
Classification  studies  using  both  short-  and  long-period 
statistics  are  in  progress. 

(2)  The  literature  search  comparing  seismic  signals  from 
nuclear  and  chemical  explosions  is  drawing  to  a  close. 
Little  useable  information  has  been  found  concerning 
teleseismic  characteristics  of  chemical  explosions.  At 
regional  distances  considerable  information  is  available. 
Chemical  explosions  in  water  appear  to  be  well  recorded 
at  teleseismic  distance.  A  report  is  being  prepared. 


(3)  The  Markov  spectral  estimate  technique  has  been 
applied  to  a  suite  of  theoretical  seismograms.  All 
the  modes,  both  normal  and  leaking,  show  up  well. 
However,  there  is  a  considerable  scatter  in  the 
estimated  dispersion.  It  is  thought  that  the  prediction 
filters  used  to  generate  the  Markov  spectral  estimates 
are  too  long.  Dispersion  estimates  using  shorter  fil¬ 
ters  are  being  run.  Preliminary  results  for  the 
theoretical  data  are  discussed  in  the  addenda  to  this 
report. 

(4)  Long-period  spectra  and  discriminants  have  been  com¬ 
puted  for  39  earthquakes  and  25  explosions  of  various 
body  wave  magnitudes.  Thresholds  are  being  deter¬ 
mined  for  usefulness  of  the  various  discriminants  for 
identifying  earthquakes  and  explosions.  The  discrimi¬ 
nants  used  are  discussed  in  the  addenda  to  this  report. 

ACTION  BY  AFTAC 

None 

FINANCIAL  STATUS 

Financial  status  for  this  contract  will  be  reported  in  the  AMSR  to 
be  submitted  the  first  week  in  July. 

Very  truly  yours, 

TEXAS  INSTRUMENTS  INCORPORATED 

Stanley  J.  Laster 

Program  Manager 
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Attachment 


2 


. 


CLASSIFICATION 


REPORT  Or  PROGRESS  AGAINST  SELECTED  MILESTONES 

_ _ (&••  ln»lrucltent  on  itvmeo  brntot*  Mltej  In) _ 

*' IT^?4?aT->  t-.  V.  fOMFONBNT 

AFTAq  Project  No.  VELA  T/0702  SEISMIC  D 

4.  NAM*  AND  LOCATION  OF  Fft*FARIN«  Activity  "  "  - - - 

lexas  Instruments  Incorporated 
Services  Group 
Dallas.  Texas  75222 


SEISMIC  DISCRIMINATION  TECHNIQUES 


».  CONTRACT  NUMBKr' 

F33657-70-C-0311 

••  NKFONT  FOR  MONTH  KNOIN* 

31  May  1970 


Code 

a 


Milestone 

b 


*  O  F  — - - 

to  Washington;  successfully  pro¬ 
cess  one  event 
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tion  to  AFTAC  for  approval 


to  use  LP  discriminants  in 
SECSP  for  data  recorded  at  LASA 


A-4  Determine  modifications  neces¬ 
sary  and  statistics  to  use  in  ap- 
plying  SECSP  to  ALPA  events 


A-5  Complete  processing  of  Eura¬ 
sian  seismic  events  (recorded 
at  LASA  and  ALPA)  using  SECSP 


A.-6  Prepare  special  report  summar¬ 
izing  results  of  SECSP  proces¬ 
sing 


Conduct  literature  search  and 
publish  report  on  differences  be¬ 
tween  seismic  signatures  of 
nuclear  and  chemical  explosions 


’Obtain  data  ensembles  with  con¬ 
tinental  paths  recorded  at  LASA 
and  A  LPA 
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ing  various  methods  of  detecting 
and  isolating  higher -mode  sur¬ 
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sing  effects  of  noise  on  various 
discrimination  statistics 
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Specify  statistics  to  use  in  eval¬ 
uating  combined  discrimination 
capability  of  LASA  and  ALPA 
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Complete  processing  oi  pre¬ 
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discrimination  capabilities 
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SECTION  I 

LONG -PERIOD  DISCRIMINATION  STUDIES 
A.  DISCRIMINANTS 

A  program  (LPSTAT)  was  written  to  calculate  the  following 
statistics  from  long-period  recordings.  (Throughout  the  following  "signal" 
refers  to  any  energy  in  the  signal  time  gate.  Similarly  for  "noise".  ) 

•  AR  parameter  for  signal  gate  (normalized  to  a 
distance  of  80°) 

•  AR  parameter  for  noise  gate  (normalized  to  a 
distance  of  80°) 

•  AR  of  signal  minus  noise 

•  Delay  time  to  the  signal  maximum 

•  Surface  wave  magnitude  estimate  computed 
from  the  signal  maximum 

•  Period  of  the  signal  maximum 

•  S/N  Ratio  (zero-to-peak  signal  maximum/RMS 
noise) 

•  S/N  Ratio  (RMS  signal/RMS  noise) 

•  RMS  parameter  for  signal  gate  (normalized  to  a 
distance  of  80°) 

•  RMS  parameter  for  noise  gate  (normalized  to  a 
distance  of  80°) 

•  RMS  of  signal  minus  noise 

•  Spectral  Ratio  for  signal  gate 

•  Frequency  at  which  the  spectral  peak  occurs 

•»  Chirp  Statistic  (amplitude  of  the  chirped  waveform 
corrected  for  distance  and  magnification) 

These  discriminants  may  be  calculated  for  all  three  components 
(vertical,  inline,  transverse).  In  addition,  the  AR  and  RMS  parameters  are 
averaged  over  the  number  of  components  input.  The  chirp  filtering  is  iterative 
where  in  the  time  duration  of  the  chirp  filter  is  successively  altered  until  the 
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maximum  crosscorrelation  of  the  filter  with  the  time  trace  is  obtained  (see 
below). 

Long -period  discriminants  have  been  calculated  from  the  verti¬ 
cal  component  of  the  Rayleigh  wave  signal  for  an  ensemble  of  39  earthquakes 
and  25  explosions  for  which  short-period  P-wave  discriminants  are  also  avail- 
able;  Two  dimensional  plots  of  these  discriminants  are  being  prepared. 

Next,  learning  patterns  will  be  formed  to  test  the  effectiveness 

of  various  combinations  of  long-  and  short-period  statistics  in  the  classification 
process. 

B.  CHIRP  FILTERING 

Several  techniques  were  investigated  wherein  the  length  of  a  chirp 
filter  is  iteratively  adjusted  to  yield  the  maximum  crosscorrelation  between  the 
filter  and  the  time  trace*  The  most  stable  technique  involves  calculating  the 
maximum  amplitude  of  the  output  for  three  different  length  filters.  The  chirp 
length  associated  with  the  maximum  of  the  function  is  determined  using  three- 
point  quadratic  interpolation.  A  chirp  filter  of  that  length  is  then  calculated 
and  applied  to  the  trace  and  the  maximum  amplitude  of  the  output  determined 
for  use  in  the  next  three  point  interpolation.  The  process  is  repeated  3  times 
and  the  maximum  amplitude  obtained  is  used  for  calculating  the  chirp  statistic. 
This  technique  has  proved  to  be  more  effective  than  using  a  chirp  length  based 
on  epicentral  distance  and  assumed  group  velocities,  but  requires  somewhat 
more  computer  time.  A  single  chirp  operation  takes  about  .  1  minutes 
(IBM  360/65)  in  the  time  domain,  while  the  iterative  procedure  takes  from  .  3 
to  .  5  minutes  (typical  chirp  and  record  length).  Both  operations  are  much 
faster  when  performed  in  the  frequency  domain. 

SECTION  H 

DETECTION  OF  LONG-PERIOD  HIGHER  MODE  ENERGY 

The  characteristics  of  the  maximum  entropy  spectrum  analysis 
technique  were  investigated  further  by  applying  it  to  seismograms  of  precisely 
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known  content.  The  purpose  of  this  was  to  perhaps  explain  some  unexpected 
results  for  the  two  events  examined  and  discussed  in  the  last  report.  The 
seismograms,  generated  under  a  previous  contract,  were  synthetsized  from  the 
solution  of  the  equations  describing  a  single  layer  on  an  infinite  half- space. 
Details  of  the  original  work  are  contained  in  the  published  report.  1 

As  the  original  seismograms  were  direct  simulations  of  an 
existing  analog  model,  scaling  of  the  physical  parameters  of  the  layer  model 
was  done  to  make  it  a  reasonable  model,  within  the  scaling  constraints,  of  the 
earth.  The  scaling  used  was  to  let  1  cm  (model)  correspond  to  20  km  (real 
earth)  and  1  microsecond  to  4/3  second.  This  results  in  the  model  below 


TABLE  1 

P  velocity  (km/sec) 

S  velocity  (km/sec) 

Density  Ratios 

L 

5.75 

3.21 

1.0 

40  km 

8.23 

4.86 

This  model  is  not  as  close  to  a  typical  continental  crust  as  one 
might  desire  because  (1)  The  crustal  velocity  is  somewhat  too  low,  (2)  The 
crustal  thickness  is  somewhat  too  small,  and  (3)  The  mantle  density  should 

be  larger  than  the  crustal  density.  Nonetheless,  the  model  is  adequate  for  the 
analysis  used. 

The  vertical  seismograms  of  four  modes:  M^,  M^,  PL^j,  and 
PL22  Were  comPuted  for  distance  from  1000  km  to  9400  km  at  400  km  incre¬ 
ments.  The  PL21  and  PL^  modes  are  the  first  and  second  order  leaking  modes 
while  the  M1  ^  and  M ^  modes  are  the  Rayleigh  and  first  higher  shear  modes. 


Laster,  S.  ,  foreman,  J. ,  Linville,  A.,  Theoretical  investigation  of  model 
for  a  layer  over  a  half-space:  Geophysics,  vol.  XXX,  No.  4,  August  1965. 
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These  four  modes  are  shown  for  a  distance  of  5000  km  in  Figure  la,  lb,  lc, 
and  Id.  Note  that  the  vertical  scale  is  different  for  each  mode.  The  composite 
seismogram,  the  sum  of  the  four  individual  modes  are  shown  in  le. 

The  composite  trace  was  convolved  by  a  filter  intended  to  simulate 
the  response  of  a  long-period  seismometer.  The  filtered  composite  trace  is 
shown  in  Figure  If  and  the  frequency  response  of  the  seismometer  is  shown  in 
Figure  2.  The  actual  LASA  L.P.  seismometers  have  peak  response  around 
0.  04  Hz;  however,  to  create  a  reasonable  seismogram  similar  in  appearance 
to  real  LASA  seismograms  the  peak  response  was  lowered  to  0.03  Hz. 

In  addition  to  the  pure  event  signal,  t.ie  composite  was  mixed 
with  random  noise,  itself  having  been  filtered  by  the  seismometer  response. 
Signal- to-noise  ratios  of  10:1  and  50:1  were  used  tsignal  zero-to-peak/noise 
rms).  These  traces  are  shown  in  Figures  lg  and  lh. 

The  seismograms  shown  in  Figure  1  are  being  individually 
analyzed  by  the  maximum  entropy  technique.  Preliminary  results  indicate  the 
following  conclusions: 

•  The  accuracy  of  maximum  entropy  technique  is 
strongly  dependent  on  the  precision  of  the  calcu¬ 
lations.  From  looking  at  the  computed  results 
it  appears  that  large  signal-to-noise  events  can 
introduce  spurious  peaks  in  the  spectrum  unless 
adequate  machine  precision  is  used. 

•  The  present  long-period  data  is  oversampled. 

The  response  of  the  LASA  seismometers  pro¬ 
duce  strong  cutoff  above  0.  1  Hz  resulting  in 
little  power  between  0.  1  Hz  and  0.5  Hz.  This 
also  causes  problems  related  to  adequate  pre¬ 
cision  and  the  introduction  of  spurious  peaks, 
when  using  the  maximum  entropy  technique. 

•  The  data  should  by  antiaiias  filtered  and  deci¬ 
mated  to  reduce  the  incidence  of  false  spectral 
peaks,.  In  addition,  whitening  of  the  data  may 
also  be  desirable.  These  techniques  are  now 
being  applied  to  the  theoretical  data  and  also 
some  real  long-period  seismograms. 


Figure  1 
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